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Abstract 

We study the stability of a vector field associated to a nearly-integrable Hamiltonian 
dynamical system to which a dissipation is added. Such a system is governed by two 
parameters, named the perturbing and dissipative parameters, and it depends on a drift 
function. Assuming that the frequency of motion satisfies some resonance assumption, we 
investigate the stability of the dynamics, and precisely the variation of the action variables 
associated to the conservative model. According to the structure of the vector field, one 
can find linear and exponential stability times, which are established under smallness con- 
ditions on the parameters. We also provide some applications to concrete examples, which 
exhibit a linear or exponential stability behavior. 
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1 Introduction 

We investigate the behavior of nearly-integrable Hamiltonian vector fields to which a dissipa- 
tive contribution is added. The vector field is ruled by two parameters, namely the perturbing 
parameter (measuring the non-integrability of the system) and the dissipative parameter (pro- 
viding the size of the dissipative term). We assume that the phase space is contracted by time 
evolution. A drift function enters the equations of motion as an unknown function; its role is 
fundamental, since it must be properly chosen in order to meet some compatibility conditions 
ensuring the existence of a normal form (compare with KAM results as in [4]). We concentrate 
on the behavior of the variables which are actions of the conservative system (i.e. setting to zero 
the dissipative parameter). We assume that the initial conditions define a resonant frequency 
for the integrable conservative system (i.e. setting to zero both the perturbing and dissipative 
parameters). Under smallness conditions on the parameters, we prove that the action variables 
stay locally bounded over a given time interval (see also [20]). The length of the time interval 
depends on the functions defining the equations of motion and, precisely, whether there appear 
also dissipative resonant terms in the original as well as in the normalized vector field. Notice 
that such a result provides a useful information concerning the transient time, namely the time 
needed to reach the attractor. 

The proof of the result is based on the construction of a suitable coordinate transformation, 
which is provided by the composition of a conservative and a dissipative change of variables. 
A similar technique, based on a non-resonant normal form, has been already implemented in 
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[8] in order to investigate a vector field of the type studied in this paper, but in the simplest case 
of non-resonant frequency. In this case, under smallness conditions on the parameters, one can 
prove that the actions stay always bounded for exponential times. As in classical perturbation 
theory, the first transformation removes the conservative perturbation to higher orders (see, 
e.g., [5]); the corresponding normal form is composed by resonant or average terms. The sec- 
ond transformation is performed to normalize the dissipative terms; the normal form equations 
defining the dissipative change of variables can be solved, provided that the drift function 
is chosen in such a way that the compatibility condition is satisfied. The final normal form 
contains just resonant and average terms up to a given order in the perturbing and dissipative 
parameters. As in classical Nekhoroshev's theorem ([20], [21], see also [2], [13]) by properly 
choosing the order of the normal form, one can determine stability bounds. The stability time 
is exponential, whenever conservative resonant terms do not appear in the equation for the 
time variation of the normalized action variable or whenever the dissipative resonant contri- 
butions are zero. In the other cases the stability time depends on the inverse of the product 
of the perturbing and dissipative parameters. The scheme of the proof, which is presented for 
a non-autonomous, time-periodic system (see also [16]) follows closely [22], where a very 
clear and enlightening proof of Nekhoroshev's theorem is given. The proof is constructive and 
it allows us to provide explicit expressions for the conservative and dissipative transformations 
(see also [3], [10], [18]). In our opinion there are several physical problems, which can be 
analyzed by our method. For example, there are many results concerning the stability of the 
(resonant) Lagrangian points in a conservative framework (see [6], [7], [14], [15], [19]), but 
none of them takes into account dissipative effects (like Solar radiation, Poynting-Roberston 
drag, Yarkowsky effect, etc.), which might significantly affect the dynamics. In this respect, 
we believe that it would be interesting to analyze these models including a dissipative effect 
by using the results contained in this paper. We provide examples of normal forms in some 
concrete one-dimensional, time-dependent model problems, which illustrate different cases 
corresponding to linear (i.e., proportional to the inverse of the product of the perturbing and 
dissipative parameters) or exponential stability times. We also provide an application of the 
theorem in order to obtain rigorous stability bounds for the previous model problems. 

The paper is organized as follows. Notations and assumptions are defined in Section 2. The 
resonant normal form Lemma and the stability Theorem are proven in Section 3. Examples of 
normal form constructions to concrete model problems is given in Section 4. An application 
of the stability theorem is provided in Section 5. 
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2 Notations and assumptions 

We introduce the ^-dimensional, time-dependent vector field, described by the equations 

X = uj{y) + ehio,y{y,x,t) + iJ,foi{y,x,t) 

y = -ehio,x{y,x,t) - ii(^goi{y,x,t) -r){y,x,tj^ , (1) 

where y G M^, {x,t) G T^+^, while the definitions and assumptions on the parameters and 
functions are the following^. 

1 . Having fixed an initial datum yo e R^, we denote by ^4 C an open neighborhood of 

2. The vector field depends on the parameters e E IR+ (perturbing parameter), fi e M+ 
(dissipative parameter); we remark that we could equally admit vector parameters, i.e. 
e e M^, jj, e Mf^, but for simplicity of exposition we present the details just for the 
scalar case e e IR+, e R+. 

3. The functions u and r] are real-analytic, f-dimensional vector functions with compo- 
nents (cjW,...,Cc;W) and (r/«, r/W). We assume that there exists a regular function 
hoo — hoo{y) such that ^^g^^^^ = <^{y)- Let ho{y, u) = hoo{y) + e R, be the unper- 
turbed Hamiltonian function associated to the conservative vector field (pL — 0) in the 
extended phase space. Let cUeiy) = i^iy)A) be the frequency vector in the extended 
phase space. Following [22] we make the hypothesis that ho is L,Af -quasi convex, 
namely there exist L, M > Q, such that for all z = {y,u) E A x ^ at least one of the 
following inequalities is satisfied: 



My)-^;|>L||^;||, ^^vv > M\\vf , V^; e R^+S (2) 



where the dot denotes the scalar product and || ■ || denotes the Euclidean norm. 

4. In the following we will use the vector field (1) in the extended phase space with t — 1 
and with u conjugated to time: 

X = uj{y) + ehio,y{y,x,t) + iJ,foi{y,x,t) 

y = -ehio^^{y,x,t) - ii(^goi{y,x,t) -r]{y,x,tj^ 

ii = -£hio^t{y,x,t) + na{y,x,t) , (3) 

where the unknown function a is introduced for later convenience (see next point). 



^The subscripts x, y, t denote derivatives with respect to x, y, t, i.e. = hy = ht = — 



dt • 
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5. We assume that /oi, goi, r] are real-analytic, £-dimensional vector functions from A x 
T'^+^ to R^, while hio, a are periodic and real-analytic from A x T^+^ to M. We re- 
mark that T], o are unknown functions, which will be properly chosen so to meet some 
compatibility requirements in order to obtain a suitable normal form (see Section 3). 

6. We assume that the vector field is dissipative and that the phase space volume is con- 
tracted by the time evolution. 

7. For a given initial datum = y{0) G A, we assume that there exists a lattice A C Z^+^, 
such that the vector function u — ui{yo) satisfies the resonance condition 

Hvo) • A; + J I = for all (k, j) E A . (4) 

We also assume that there exists K E Z^, a > and a subset D C A, such that for any 
y E D the. following condition is satisfied: 

Hy) ■k + j\>a for all (k,j) E Z^+^\A , + \j\ <K , (5) 

where for k = {ki, ke) E 1} we define the norm \k\ = \ki\ + ... + \ke\. 

8. We refer to 77 = ri{y, x, t) as the drift vector function with components {rj'^^^y, x, t), 
j]^^'' {y, X, t)) that we expand as 

00 m 

V^'\y,x,t) = EE^S-.-+i(^'^'^)^'/^'""' ' k = l,...,i. 

m=0 j=0 

In a similar way we expand a in (3) as 

00 m 
m=0 j=Q 



Remark 1. We remark that for = the equations (3) reduce to the conservative vector field, 
associated to the nearly-integrable Hamiltonian function in the extended phase space 

H{y, X, u, t) ^ hoo{y) +u + ehwiv, x, t) , (6) 

where cu{y) = dhoo{y)/dy. Notice that the Hamiltonian (6) is integrable as far as the perturb- 
ing parameter is zero, i.e. e = 0. Since the vector field (3) is dissipative (see assumption 6.), 
the energy associated to (6) is decreasing with time. 

We adopt the following notations and definitions for functions and norms. 
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i) Integer subscripts denote the order in the perturbing and dissipative parameters, i.e. Fij 

denotes a function of order £*//^. 

ii) For a function f — f{y,x,t) and for any positive integer K, we decompose / as 

/(y, X, t) = f{y) + /("'■'^^Hz/, ^, t) + f'''-''\y, X, t) + /(>^)(y, t) , 

being, respectively, the average, the sum over the non-resonant components with Fourier 
modes less or equal than K, the projection over the resonant space defined by the lat- 
tice A excluding the origin with Fourier modes less or equal than K, the sum over the 
components with Fourier modes greater than K, namely: 

(fc,i)€A\{0}, |fc| + b|<i^ 



where i — v — 1 and fkj are the Fourier coefficients. 

Hi) We say that a function is of order k in e and /x, in symbols Ok{£, A*), if its Taylor series 
expansion in e, ji contains powers of e'^ji^ with i + j > /c, i > 0, j > 0. 

iv) We denote by Cr^iA) the complex neighborhood of A of radius ro, namely 

Cr^A) = {y e : \\y - y^H < ro for all yA^A}. 

Moreover, let Cso(T^+^) be the complex strip of radius sq around T^+^, namely 

Cs,{T'+^) = {{x,t) e C'+^ : max \Q{xj)\ < so , m)\ < so} , 

i<j<e 

where denotes the imaginary part. 

v) Denoting the Fourier expansion of a function / = f{y, x, t) as 

/(t/,x,t)= Y hjiyy^''^^''^ ^ 

(k,j) 
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we introduce the norm 



sup Yl 

For a function g = g{y) we define \\g\\ra = ^'^VyeCroiA) \\9iy)h where || • || denotes the 
Euclidean norm. For a vector function / = (/i, .., /^) we define 



^•0,50 



3 Bounds on the variation of the action variables 

In order to bound the variation of the action variables, we implement a change of coordinates 
such that the vector field (3) is transformed to a resonant normal form, up to a suitable order 
A^. To this end we introduce a change of coordinates close to the identity and leaving time 
unaltered, say 

(y,X,C/,f)=sW(y,x,ii,i), {Y,U)eR'+\ {X,t)eT'+\ (7) 

where S'^^) depends parametrically also on e, ji, E'^^^ = E'^^^ [y, x, u, t; e, jS) with S'^^) [y, x, u,t;0, 0) 
Id. Let K be as in (5); in the forthcoming Resonant Normal Form Lemma we aim to determine 
the transformation of coordinates (7), so that (3) takes a resonant normal form of order N, that 
we write as 

X = nf\Y) + F^'^^''\Y,X,t) + FM+iiY,X,t) + F^>''\Y,X,t) 

Y = G^'''^''\Y,X,t) + GN+i{Y,X,t) + G^>''\Y,X,t) 

jj = H^-'^''\Y,X,t) + HN+,{Y,X,t) + H>''\Y,X,t), (8) 

where Q^J^^ : ^ is the normalized frequency, related to oj{Y) by 

AT N N-j 

i=l j=l i=0 

where Qij{Y) are known vector functions; F^'^'-^\ G^'^'-^\ j{{r,<K) ^^ve Fourier components 
belonging to the resonant lattice A\{0} with F^'^'-^^ depending on both e, fj,, while G^'^'-^\ 
}{ir,<K) (iepend only on e; F^^i, G^^i, H^+i are vector functions of order On+i{s, jj); 
pi>K)^ Qi>K)^ jji>K) (denote functions with Fourier modes greater than K. 
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Similarly to [8] we decompose the coordinate transformation S^^^ as the composition of two 
transformations S^^^ (conservative part) and S^^^ (dissipative part): 

(Y,X,U,t)^EfKEi^\y,x,u,t). (9) 

Setting (y, x, u, t) = E!i'\y^ x, u, t), the conservative transformation si^^ is defined through 
a sequence of generating functions close to the identity, say tpjo = ipjo{y, x, t), j = 1, N, 
such that 

N 

X = x + ^i^jo,y{y,x,t)e^ = x + 'ipl^\y,x,t) 
i=i 

N 

y = y + ^'^jo,x(y^^^'t)£^ = y + 'ijj^^\y,x,t) 

AT 

u ^ u + ^ iljjo,t{y, X, t)e^ = M + il^f'^ {y, x, t) . (10) 
j=i 

Notice that we can assume that the functions V'jo (as well as ajk, i3jk, ^jk in (12) below) do not 
depend on u, since the functions appearing in (1) (or equivalently in (3)) do not depend on u. 
We denote the inversion of (10) as 

X = x{y,x,t)^x + r^'''^\y,x,t) 
y = y{y,x,t)^y + T^y'''\y,x,t) 

u = ii(y,x,i)={i + r("'^)(y,x,i). (11) 

The dissipative transformation S^^^ is defined by introducing suitable functions with zero 
average over x and t, say a^^\ /3'^^\ 7^^^ defined through series by the coefficients aji, /3ji, 
Iji, 3, i e such that 

AT i 

X = x + ^^aj^i-j{y,x,t)e^ii'~^ = x + a^^\y,x,t) 

N i 

N i 

i=o 3=0 

with Q;io(y, X, t) — /3io{y, x, t) = 7jo(^, x,t) = for any i > 0. An iterative explicit construc- 
tion of the vector functions ipjo, aji, /3ji, 7^^ will be given within the proof of the Resonant 
Normal Form Lemma stated below. 
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Remark 2. The normal form equation defining the generating function ipjoiy, 3:,t) at order j 
is given by 

i^{y) i^jo,x{y^ + ^jo,t{y^ ^) + LfQ'-^\y, X, = , 

for a suitable known function L^^^'-^^ (y, x, t) with zero average over (x, t) and not containing 
resonant terms, say 

{fcj)ez^+i\A, |fc|+b-|<i^ 

This equation can be solved provided lu = u}{y) satisfies a non-resonance condition of the 
form 

uj{y) -k + j^O for all (k,j) e Z^+^\A , + |j| < K , 

which is guaranteed by (5), provided e satisfies a smallness condition. Analogously, the dissi- 
pative normal form provides an explicit construction of the functions a'^^\ /3^^\ 7*^^^, thanks 
to a suitable choice of the drifts 77, a and to the assumption (5). More precisely, once expressed 
in terms of the new variables (F, X, t), the functions Pji must satisfy a normal form equation 
of the form 

uj{Y)P^UY, X, t)+(3ji,t{Y, X, t)+N^r-''\Y, X, t)+Nji{Y)+N^l'^''\Y, X, t)+riji{Y, X,t)=0, 

(13) 

for some known function Nji = Nji + jvjf^'-^-' + Njl'-^^ + Njf^^; therefore, equation (13) 
can be solved provided the drift components tjji^Y, X, t) are chosen as the opposite of the sum 
of the average and of the resonant parts: 

ri,,{Y, X,t)^- (^Nj,{Y) + N^I'^''\y, X, t)) . 

An analogous relation holds for and aji. This explains why the drift must be properly 
defined in order to be able to build the coordinate transformation (7). This is not unusual, but it 
happens also in KAM theory (see e.g. [4]). We proceed now to state the Resonant Normal Form 
Lemma, which extends the Normal Form Lemma of [8] to the resonant case of a frequency 
vector satisfying (4), (5). 

Resonant Normal Form Lemma. Consider the vector field (1) analytic in the complex ex- 
tension Cr„{^) ^ C's(,(T^^^) for some tq, sq > 0. Consider the extended vector field (3) on 
AxRx T^+i. For a given lattice A C Z^+\ let yo E A, K e Z+, D C A, a > be such 
that (4) and (5) are satisfied. There exist suitable drift functions r] = ri{y, x, t), a = a{y, x, t), 
and there exist Eq, /Iq > depending on tq, sq, K, a and on the norms of u, h, f, g, such 
that for any {e, /i) < {eo, /io), one can construct a change of variables close to the identity, say 
S(^) : >1 X R X T^+i ^ M^+^ x T^+^ with (Y, X, U, t) = S^^) (y, x, u, t), being (Y, U) e M^+\ 



9 



{X, t) E T ^+1^ ]V G Z+, which transforms (3) into a normal form of order N as in (8). Let 
Rq < tq. So < So; having set A = max(£, //), the normalized frequency is bounded by 

\\^T^-uj\\n,<C^X, (14) 

where C^j is a positive constant depending on ro, N and on the norms ofu, h, f, g. Denoting 
by Uy the projection on the y -coordinate, one gets 

||n,(sf)osf))-/4<c,A, (15) 

for some positive constant Cp depending on ro, sq, and on the norms of cu, h, f, g. With 
reference to the normal form (8), one has the following estimate 

||G^'''-''^lk,5o + \\Gn+iU,So + l|G(>^)||iJo,5o < + , (16) 

for some constant Cy and having bounded || G*^^'-^) \\ro,So by A G, where G and Cy depend on 
ro, So, N, K and on the norms of u, h, f, g. Choosing^ N — [Kto/\ log A|] for some tq > 0, 
one obtains that (16) becomes 

IIG^'-'^^^lkA + \\Gn+iWso + l|G(>^^l|iio,5o < AG + CyAe-^-° . (17) 



Before giving the proof of the Lemma, we provide the statement of the main result, namely 
a bound on the variation of the variables which are actions of the conservative system. The 
following Theorem will be obtained through the Resonant Normal Form Lemma under the 
resonance condition (4) and the quasi-convexity assumption (2). Let us write the normal form 
equations (8) using the following notation: 

X = u;{Y)+epi^''\Y,X,t)+fis^^''\Y,X,t) + FN+i{Y,X,t) + F^>''\Y,X,t) 
Y = -ep^^''\Y,X,t) + GN+^iY,X,t) + G>^\Y,X,t) 

U = -epi-''\Y,X,t) + HN+i{Y,X,t) + H>''\Y,X,t), (18) 

where p^x^K P^y^\ pi~^^ (independent of fi) are the resonant contributions stemming just 
from the conservative transformation, while s^-^^ (depending on and e) represents the reso- 
nant part coming from the dissipative transformation. 

Theorem Consider the vector held (1) defined on A x T^+^, satisfying the quasi-convexity 
assumption (2). Letyo E A, K E Z+, D C A, a > be such that (4) and (5) are satisfied. 
Assume there exists Eq, i^o, such that for (e, /i) < (^o, l^o), the Resonant Normal Form Lemma 
holds. Let tq, Cp, A, ro, sq as in the Resonant Normal Form Lemma. With reference to (18), 
we have that: 

^The choice of N is motivated as follows. The relation = e~^^° implies A'^logA = —Ktq, namely 
N = [Kto/\ log A|], where [•] denotes the integer part. 
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i) ifp^i ' = or .s(^^^) = 0, then there exist pi > 0, Cq > 0, such that \\y{t) - |/(0)|| < 
2Cp\ + pi for t < Ti = Coe^'^", where Cq depends on M, tq, sq, K, N and on the 
norms of uj, h, f, g, while pi depends on the above and on m, X, A; 

a) ifp^^^^ and s^^^^ 0, then there exist p2 > 0, > 0, C'^ > 0, such that if 
t<T2 = min(C^e^^", then \\y{t) - y{0)\\ < 2CpX + p2, where C'^, C^' depend on 
M, To, So, K, N and on the norms of to, h, f, g, while p2 depends on the above and on 
m, X, A. 

Remark 3. The above theorem is stated in terms of the functions p^x^^ ™^ s^-^^ appearing 
in the normal form equations (18); in order to decide which of the conditions i) or ii) of the 
Theorem is satisfied, one needs to know the explicit expression of the functions /oi, goi, ^lo 
appearing in the vector field (1), tracing the resonant terms which could generate p^x^^ ^^'^ 
si<K) by means of an explicit construction of the normal form or by means of a tree algorithm 
(see, e.g., [9], [11], [17] and references therein). 

Remark 4. The Theorem states that in the non-resonant case (compare with [8]), as well as 
whenever the dissipative contribution to the resonant normal form is zero (at least up to the 
normalization order), one finds a variation of the actions on exponential times; otherwise, there 
appears a fast drift of the actions on linear (in e p) times. 

Proof of the Resonant Normal Form Lemma. By induction on the normalization order we 
prove that we can construct a normal form of type (8) by means of suitable transformations 
as in (10) and (12). First we prove the statement by constructing the first order normal form 
using the conservative and then the dissipative transformation; next, we proceed to construct 
the conservative and dissipative transformations at the order N. For sake of clarity, we split 
the proof into four separate steps, referring, respectively, to the first order conservative and 
dissipative normal forms, and to the iV-th order conservative and dissipative transformations. 
Since the conservative transformation is standard, we omit some details. 

Step 1: Conservative transformation for = 1. 

We start by implementing the first order transformation 

y = y + ei)iQ^^{y,x,t) 
u = u + eipio,tiy,x,t) , 

where V'lo = ipioiv, x, t) is an unknown function. Let fo < ro, < sq, sq 
can invert (19) as 

X = x + sr^''''^\y,x,t) 



(19) 

= So — 5o; then we 
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u = u + er^'''^\y,x,t) , 

for suitable functions F^^'^^, T^V'^^ and F^"'^), provided the following smallness condition on 
the parameters is satisfied (compare with Appendix A): 

70 ||^io,y||fo,.o e''°^o"'^ < 1 • (20) 

Using (19) and (3), we obtain that the conservative normal form is achieved whenever one can 
determine V'io(y, x, t) such that 

t^2/(y)V'io,x(y, X, t) + uj{y)ilJio,yx{y, x, t) + ^w,yt{y, x, t) + h!{Qf^\y, x,t) = 
t^(^)^/'io,xx(^,5,t) + ^/'io,xt(^,5,t) + Mo^^^l^,^,^) = 
'^(^)^io,te(^,5,t) +^io,tt(y,5,t) + /iSo|t'"'^Hy,5,^) = 0.(21) 

Let us define 

Equations (21) are equivalent to take the derivatives with respect to y, x and t of 

oj{y)ilJio,xiy, X, t) + V'io,t(y, X, t) + h^w'~^\y^ x,t)^0 . 

Expanding tpio and h^^^'-^^ into Fourier series, one obtains that ipio is given by the expression 
(independent of u): 

^,o(y,x,t)=z ,l^. e-(^-^+^-*^ (22) 

{k,j)ei.(+^\A,\k\+\j\<K 

This function is well defined, since the zero and small divisors are controlled as follows. The 
second of (19) can be inverted as y = y + £it!^^'^)(y, for a suitable function R^^'^^ — 
R^^'^^ (y, X, t), provided that for fg < Tq one has (see Appendix A) 

70s\\i'io,Aro,so^:r^ <l , (23) 
being \\f'Q,so ^ ||^io,a;||fo,so- Then, the divisors appearing in (22) are bounded by 

\uj{y)-k + j\ >a-eK\\R^y''^\y^,J\ujy\\ro > | , (24) 



provided 



£ < ^r.nr^U.UU r^-ir- ■ (25) 



2K\\R(y'^)y^ 



soU^yWro 
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Step 2: Dissipative transformation for iV = 1. 

We define the first-order dissipative transformation as 

X = x + Q;oi(y,x,i)// 
Y = y + poi{y,x,t)iJ, 

U = u + -foi{y,x,t)ii , (26) 
for unknown functions ctoi, Poi and 701. Let us start by inverting (26) as 

X = X-Q;oi(y,X,i)// + 02(/i) + A("'i)(y,X,i)/i 

y = y-^oi(i^,^,i)/^ + 02(//) = y + A(^'i)(y,x,i)// 

It = U-^oi{Y,X,t)pL + 02{pL)^U + A^''''\Y,X,t)pi (27) 

for suitable functions A^^'^^, A^^' i) and A^"'^) provided the following smallness conditions on 
the parameters are satisfied (see Appendix A): 

70 ||Q;oi||fo,So < 1 

70 (||^Ol||fo,so + l|/^01,a;||fo,so l|Q;oi||fo,5o) 5~ A* < 1 

70 (||7oi||fo,so + ll7oi,x||fo,so ||«Ol||fo,so) JT < ^ ' (^^) 

where = 5o/2, Rq < fo and being ||A(^'^)||-g -^.^^ < ||aoi||fo,So- Through (26) and (27) we 
can express X, y as a function of X, Y; the normal form is obtained assuming that aoi, /3oi 
and rjoi satisfy the following equations: 

iu{Y)aoiAy, X, t) + aoi,t{y, X, t) - iJy{Y)/3oi{Y, X, t) + ft'-''\y, X,t)^0 
u;{Y)l3oi,.{Y,X,t) + l3oi,t{Y, X,t) 

- giT-'^hY, X, t) - M{Y) - gi'{-''\Y, X, t) + voi{Y, X,t)^0 
uj{YhoiAY,X,t) + 7oM(F,X,t) + aoi(y,X,t) = 0. (29) 

Since aij, Pij, '^ij have zero average and they do not contain resonant terms, the system of 
equations (29) can be solved, provided that we choose r]Qi{Y, X, t) as 

r/oi(y, X, t) = -goiiY) + gt{^''\Y, X, t) 

and that we set aoi = as well as 701 = 0. Setting fi^^-* = (jj{Y) + ehio,y{Y) + /x/oi(y), the 
first order normal form can be written as 

X = ^lf+eht,f\Y,X,t)+^ft^^\Y,X,t) 

+ eh^^^fiY, X, t) + X, t) + F2(F, X, t) 

Y = -eht^J'\Y,X,t)-eh^>^\Y,X,t)-^,g^^,''\Y,X,t) + G,{Y,X,t) 

U = -eht^^''\Y,X,t)-eht^f{Y,X,t) + H2{Y,X,t), 
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where F2, G2 are functions of order 02{e, ji) and H2 is a function of order 02(e). We remark 
that the solution of (29) involves small divisors of the form uj{Y) ■ k + j with {k,j) G Z^"*"^ 
and |A;| + |j| < K. Using the same argument as in (24), the small divisors are bounded by a/4 
provided that the following smallness condition holds (compare with Appendix A): 

< .r..n II 11 ■ (30) 

Step 3: Conservative transformation for the order N. 

Assume that the Lemma holds to the order N—1. We introduce the conservative transformation 
to the order N as 

N-l 



N-l 

N-l 

u = ^^jo,t(y,a;,t)£^ +^Aro,t(y,a;,t)£^ = M + ^J^-'(y,a;,t) , (31) 



where for 1 < j' < — 1 the functions ipjo{y., x, t) are assumed to be known. We can invert 
(31) as 

X = x{y,x,t) 

N 

^ X + ^foiy, t)e^ - ^mM t)^"" + ON+i{e) = x + r^-'^\y, x, t) 

y = y(y,x,t) 

N 

^ y + J2 (y^ S:, ty + ^ , t)e'' + 0^+1 (e) = y + r(2''^) (y, t) 

i=i 

u = u{y,x,t) 

N 

^ u + Y, (y, X, t)e^ + i^mM ^> + ON+i{e) = u + r("'^)(y, i) , (32) 

provided that (see Appendix A), choosing fo < ro, 5o < so: 

70||^ri|.o,.oe'^°^o"' < 1, (33) 
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)>K 



being F^q^ F^q^ F^q\ 1 < j < A^, known functions. We proceed to compute x, y, u as a 
function of x, y, t and after expanding in Taylor series, we obtain 

'N+1 

y = x,t) + i, {m-iAy, ^, t)^""'' + - + voAy, ^, ^)/^'^"') + oj 

where O^^^ is a compact notation to denote terms of order On+i{£,I^) and/or containing 
Fourier components greater than K; the functions F^^'-^'-^\ G^^'-^'^\ }j{^,<k,n) known, 
contain Fourier components up to the order K, contain orders in e and /x up to the order and 
they are at most linear in fi. Using (31), (34) and the inductive hypothesis, the conservative 
normal form at the order N is obtained once the function ip^o satisfies the equations 

^v{y)'^m,x{y^ ^) + ^(y)'^No,yx(y, x, t) + 'il^m,yt(y, X, t) + F^^^'''-^'-^\y, x,t) = 
i^(y)'^m,xx(y, X, t) + il^mMy^ - G%Q'''-^'-^\y, x,t) = 
u{y)i^m,tx{y, X, t) + i^m,tt{y, X, t) - H^^^'''-^'-'^\y, x,t) = , 

(35) 

where F^q^^'-^'-^\ G%q^'-^'-^\ H^^^^'-^'-^^ are the non-resonant parts of known func- 
tions F(2.<^.<^), G^'^,<K,<N) H{2,<K,<N) that we decompose as 

AT AT 

F(^'^^'^^\y, x,t) ^ Yl Fj}-''\y)e^ + E Fj^'-'^'-'^'iy, x, t)e^ 

j=i j=i 

N N-l 

+ E ^o"-"'-"^(^, t)e^ + ^ E F^l'-"'-"\y^ t)s^ (36) 

j=l j=0 

(and similar for the remaining functions). From the Hamiltonian structure it can be easily 
recognized that F^^'^'-^'-^\ G%q""-^"-^\ H^^q""-^"-^^ are, respectively, the derivatives 
with respect to y, x, t of the same function, so that equations (35) uniquely define the solution 
ipNoiy, x, t). We are finally led to the following conservative normal form: 



N N-l 

j^{2,<K,<N)i~ ~ ^^^j ^ ^ 

'N+1 



i - ^f\me) + Y.Fir^-''^-''\mx,^^^^ 

3=1 j=0 

N N-l 

~ = E^fo"'"''"^(^'-'^)^^+/^E^fi''"''"^(^'^ 

j=l 3=0 
N N-l 

~ = Y: H^r'^^'^'^'iy, X, t)e^ + n E Hri-'''-''\y, + 0>nIi , (37) 

j=l j=0 
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where 

N 

which implies that — ^^^H < CcS for a suitable constant Cc- The normal form equations 

can be solved, provided that the small divisors taking the expression u!{y) ■ k + j, for k G Z^, 
J G Z with I A; I + \j\ < K, are controlled by a non-resonance condition, which is guaranteed 
whenever (see Appendix A) 

a 

' - 2X||i?(^.^)||,,,,J|c.,||.„ ' ^^^^ 

where R^y'^^ is the function inverting the transformation, namely y = y + eR^^'^^ (y, x, t), and 
fg < fo. The inversion can be performed provided (see Appendix A) 

70||V'f^l|.-o,.o-^ < 1 (39) 



with||i?(^'^)||,.,,„ < llV^i^ 



Step 4: Dissipative transformation for the order N. 

We consider the transformation (12) at the order A^, which can be inverted as 

X = x{Y,X,t) 

N-2N-l-i N-l 

- X + Y. Yl a,,{Y,X,t)e'^^-J2'^^,N-^iy,X,t)e'^l''-' + 0r,+l{e,^l) 

i=0 j=l i=0 

y = mx,t) 

N-2N-l-i N-l 

^ y + J2Yl %(>^,^,^)^V-^A,iv-^(F,x,^)£>^-' + o^+l(£,/.) 

i=0 j=i 1=0 

u = u{Y,X,t) 

N-2N-l-i N-l 

- U+Y Yl ^^^■(^' ^' - Y ^i,N-i{Y, X, i)£>^-^ + ON+i{e, fi) , (40) 

i=0 j=l i=0 

for suitable known functions aij{Y, X, t), bij{Y, X, t), Cij{Y, X, t), provided that the parame- 
ters satisfy (see Appendix A): 

70(||/3(^)||.o,.-o + ll/3f^l|.-o,.-oll«^^^l|.-o,.-o) < 1 

ro — rto 

70(|UW||- - + lUWll- - lUWll- - ) I < 1 (41) 

'"-"Vll' \\ro,so ^ W Ix lko,soll" \\ro,so) ~ r> ^ ) v^^^ 

'^O ~ -n-o 
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where 60 = so/2, Rq < fo. In order to determine the unknown functions ao,N, 0:^-1,1, 
/3o,iv, Pn-1,1, 1o,n, ■■; 7n-i,i, Vn-1,0, Vo,n-i,^(^n-i,o, (^o,n-i, using (37) and (40) we 
express x, y in terms of X, Y and we compute X, Y using (12), (37), (40) as 

/N-l \ N-1 

X = U{Y) - Uy{Y) f^^,N-^{Y, X, + u{Y) <^^,N-^Ay, X, t)6'fi' 



.N-i 

, i=0 / i=0 



iV-1 



=0 



Af-l 

Y = yi(r]N-i,i{Y^ X, t)e''-' + ... + rjoMY, X, i)//^"^) + u{Y) ^ ^,N-iAY^ X, i)£>^-^ 

i=0 

AT-l 

+ J] A,iv-i,t(i^, X, i)£>^- + G'(3'-'^^'^)(y, X, t) + G(=^'^)(y) 

i=0 

AT 

1=1 

TV-l 

tl = /x((7Ar-i,i(>^, ^, tje"^"' + - + <^o,N{y, X, i)Ai^-^) + u{Y) J] %,Ar-i,.(r, X, i)£>^-* 

i=0 

TV-l 

+ E 7i,iv-i,t(l^, X, + //(3.-.<^.^)(y, X, i) + ^(^■^)(y) 

i=0 
TV 

1=1 

where F^^'^^'-^'^\ G(3."'-.<^.^)^ fj{3,nr,<K,N) denote known non-resonant functions that we 
can expand as 

N-l 

piS,nr,<K,N)^Y,X,t) = ^ '"^^^^^ ^> 

i=0 

and similarly for G^3,nr,<K,N) ^ff{3,nr,^K,N). p{3,r,<K,<N)^ G^3,r,<K,N)^ jj{3,r,<K,N) ^^^^^^ ^^^^^ 

resonant functions; F^^'--^\ G^^'-^\ ^(3,^) denote the average terms. Recall that due to the in- 
ductive hypothesis, the functions aij, Pij, jij, rjij, Uij, determine a normal form up to the order 
with 0<i+j<A^ — 1. The normal form at the order N is obtained by imposing that 

otij, Pij, %j, Vij^ ^ij satisfy the normal form equations 

-ujy{Y)l3i,N_i{Y, X, t) + uj{Y)ai,N-iAY, X, t) + a^^N-iA^, ^> ^) 
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u;(y)A,Ar-i,.(y, X, t) + Pi,N-iAy, X, t) + cf^rf x, t) + g[%'^/''^''\y, x, t) 

+G!yi(i^) + r/,Ar_.(r,x,t) = o 

uj{YH,^_,,AY, X, t) + ^,,N-^,t{y. X. t) + Hl%-:f''^''\Y, X, t) + X, t) 

+H^,{Y) + a,,r,-i{Y, X,t) = Q (42) 

for < i < — 1. Equations (42) can be solved provided uj{Y) satisfies a non-resonance 
condition, which is guaranteed by (see Appendix A) 

4Ala;,||,J|/3(^)||,„,,„ <a, (43) 

where we intend that Y = y + /3^^^ {y, x, t; e, ji). From the second and third of (42), we get 

r/,^_,(F, X, t) ^ G^i^l^) + Gf0^^^\Y, X, 

X, t) ^ Hj;f_\{Y) + Hj;0^'^\Y, X, t)e^^^-^ . 

Setting 

N N N-j 

flf\Y; s, /.) ^ ij{Y) + J2 ^io{Y)e' + ^^.■(^)^V , (44) 

i=l j=l i=0 

with Qio = fIq'-^\y) and fli,N-i = the normal form is finally given by 

N N N-j 

X = nf\Y;e,,) + J2F^t''''''''\Y,X,t)e^ + J2J2 

i=l j=l i=Q 

+ FN+,{Y,X,t) + F^>''\Y,X,t) 

N 

y = 5]Gg'''-'''-'^^(>^,^,^y+Giv+i(i^,x,i)+G(>^)(y,x,i) 

i=l 

N 

i=l 

where F^v+i, G^v+i are ON+i{e,^i), Hn+i is order ON+i{e) and F^>^\ G>^\ H^>^^ con- 
tain only terms with Fourier index greater than K. The normal form (8) is recovered with an 
obvious identification of the functions F^^'-^\ G^'^'-^\ j{{r,<K)^ -pj^g smallness requirements 
on e, II, say e < Eq, /i < /iq, are needed to guarantee the non-resonance condition (see (25), 
(30), (38), (43)) and the inversion of the transformations (see (20), (23), (28), (33), (39), (41)). 
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The estimate (14) holds true, due to the definition of il^^^ in (44). The estimate (15) follows 
from the fact that (9) is close to the identity up to first order. 

Due to the exponential decay of the Fourier coefficients (compare with Lemma B.l of Ap- 
pendix B), we can bound for some tq > as 

||G(>'')|Uo,5o<C'GAe-^^ (46) 
for a suitable constant Cq- On the other hand we can bound Gn+i in (45) as 

||G'^+i||ro,5o<CgA^+\ (47) 
for a suitable constant Cq- Finally, from the second of (45) we obtain: 

N 

II y^,Gf/'-^'-^\Y, X,t)e^\\R^^So + ll<^iV+l||Ro,5o + IIG'^^'^^IUcSo 
i=l 



Kto 



having defined AG as an upper bound of sup^<^|j || Xlili G^^'^'-^'-^\y, X, t)e'^\\Rg^So- Choos- 
ing = [Ktq/\ logA|], we obtain (16) and (17) with Gy = Gq + Cq- This concludes the 
proof of the Lemma. □ 



Proof of the theorem. The distance between y(t) and the initial condition y(0) for i > can 
be bounded by the sum of the following terms: 

Mt) - 1/(0)11 < \\y{t) - Y{t)\\ + \\Y{t) - F(0)|| + ||F(0) - |/(0)|| . (48) 

By the estimate (15) of the Resonant Normal Form Lemma, one obtains 

hit) - Y{t)\\ < A , ||y(0) - y(0)|| <G,X. 

By the second of (8) and by (16), one gets 

\\y{t) - ym < £ {\\g^'''-''^\u,so + \\Gn+i\u,so + \\G^^''^\u,so)ds 

< XGt + GyX^+H , 

which indicates that the action variation takes place on linear time scales due to the term 
XGt, while exponential times are associated to the term Cy A^+^ t. We remark that G — 
corresponds to the absence of resonant terms in the normal form for Y, Notice that the case of 
non-resonant stability estimates given in [8] is recovered whenever also the resonant terms in 
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the X variable are zero. Let us start with the case G — 0; for a suitable pi > that we write 
as pi = Cp\ for some Cp > 0, let 

t < ^ e^-° . 

Finally, setting po = (2Cp + Cp) A, we obtain the following variation of the original action 
variables on exponential times: 

\\y{t)~ym\<Po fort<Coe^^°, 

having defined Co = Cp/Cy- This result is in agreement with statement i), once IS 
identified with —ep~^^^ (compare with (18)). 

Next, we study the case G 7^ 0; to this end, we compute the variation of the energy (i.e. the 
Lyapunov function, see e.g. [1]), which we intend to be defined as the energy function which 
is preserved whenever p = 0. Let us write the normal form equations (45) using the following 
compact notation as in (18): 

+ F(>^)(r, X, t- £, p) + X, t- £, p) 

Y = -£pi-'^)(y,X,i;£) + G'(>^)(y,X,t;£,p) + G'^+i(y,X,i;£,p) 
U = -ep?''\Y, X, t; e) + H>''\Y, X, t; e, p) + Hm+i{Y, X, t; e, p) , (49) 

where we have indicated also the dependence on the parameters and we have identified the 
functions as follows: 

N N 



ep'y^''\Y,X,t;e) ^ J2^,o{Y)e^ + J2F^t'-'''-"'\y,XM 

i=l i=l 
N 



i=l 
N 

spP)(y,x,t;e) ^ -j;i/r'^-'^-)(y,x,t)s' 

i=l 

N N-j N N-j 

^s^^^\Y,X,t;s,^) ^ ^^i^f'^^'^^n>^,^,t)^V + EE^^^(^)^V- 

j=l 1=0 j=l i=0 

For p = equations (49) reduce to 

X = uj{Y)+ep^^''\Y,X,t;e) + F^>''\Y,X,t-e,0) + Fj^+i(Y,X,t;e,0) 

Y = -ep^}''\Y,X,t;e) + G^>''\Y,X,t;e,0) + GN+i{Y,X,t;e,0) 

if = -ep'i^''\Y,X,t;e) + H>''\Y,X,t;e,0) + Hj,+,{Y,X,t;e,0). (50) 
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Due to the Hamiltonian character of the equations of motion for fj, — 0, there exist vector 
functions A^^^\ B^+i, such that 

j^i^K) ^ F{>K)^Y,X,t;e,0) , A^^""^ = -G>''\Y, X,t- e,0) , A^-^^ - -H>''\Y, X,t- 6,0) 
Bn+\,y — Fn+i(X^ ^1 ^) 0) 5 Bn^i^x = —GN+iiy, X, t; s, 0) , -Bat+i,* = —Hp^+i(Y, X, t; s, 0) , 

so that we can recognize (50) as Hamihon's equations associated to the following Hamiltonian 
function in the extended phase space with i = 1: 

H{Y, X, U, t) = hoo{Y) + U + £p(^^)(y, X, t) + A(>^)(y, X, t) + Bn+i{Y, X, t) , 

where /iqo is such that ^^gy^^ — cu{Y). Let us fix the energy level H — E for some real 
constant E; taking into account the complete equations (49), we obtain that the variation of E 
for // 7^ is given by (for simplicity we omit the arguments): 

^^e^,p^i''^s^^^^ + C^^^ + D>^\ (51) 
at 

with 

Cn+i = u(Y)Gn+i + epi^^^GN+i + ep-^^^FN+i 

+ (4>^^ + Bn+^^x)Fn+i + Bx+iA^ + epf''^ + + F(>^)) + i/jv+i + ^iv+i.t 

where now the functions Fjv+i, G^v+i, Hn+i, F>^\ G^>^\ H>^') depend on (F, X, t; e, fx). 
Denoting by AE = E(t) - E(0), we obtain 

where A/iqo + AU = hoo{Y{t)) - hoo{Y{0)) + U(t) - U{0) and similarly for the other 
quantities. Recalling (51) and setting ho = hoo + U, we get 



\Aho\ < |AE|+£||Ap(s^^)|U„,5o + ||A^(>^)|Uo,5o + ||A5Ar+i|ko, 



So 



where 

,dE 



A^l <\^\t< {si,\\p^i''^\\j^,sM-''^\U,So + \\Cn+i\U,So + l|i^^>^^lko,5o) t 
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sup II II ^ ^ ■ 



We denote by m an upper bound on the Hessian of /loo(y) and let m be an upper bound of the 
Hessian in the normalized variables, which we can define as m = m + || ^^||ro ||-D*^^'^^ llr-o.so' 
having expressed the link between new and old variables asY — y + D^^^^^ [y, x, t). Then, we 
have: 

Assume that the frequency uj(.{y) = {u}{y),l) is close to exact A-resonances (compare with 
[22]) by a quantity 5 > 0, namely if i?A = {Q e R^+^ : Q • n = for all n e A}, then 
min^^gR^ lli^ely) — ^11 < S. Setting Z = {Y,U), assume that \Z{t) — Z{0)\ < r for some 
r > with 6 + mr < Rq. Let IIa be the orthogonal projection on A; by the mean value 
theorem we obtain 

\uje ■ HaAZI < ||nAa;e||j?o ||PAZ|| < {5 + mr) \\AZ\\ 

\Ue-iId.-UA)AZ\ < t\\uJe\\Ro (||G'(>^^ ||fl„,5o + ll^iv+l |k,5o + ll^^^""^ lk,5o + II^^AT+l ||ilo,^^ 

Moreover 



Jo 



Aho = uje-AZ+ I (1-s) '.^-^AZ ■ AZ ds , 

so that, in the region where the M-convexity (2) holds, one has 

^||AZ||^ < \uje ■ AZ\ + \Aho\ , 

where (similarly to m) we can set M = M — || ||ro H-D*-^'^-* (y, ||ro,so' that one has 

^?|^^-^>M||^||2, V^;eM^+^ 

Finally, we have 

^\\AZf < {6 + mr)\\AZ\\+t\\uj,\\R,,(^\\G^>''^\\R,,So + IIG^+ilji^c^o 

+ ||-H'^^^^||i?o,5o + ||-f^iV+l||iio,So) + ^fJ'\\P^X^^\\Ro,So ||s^-^^||ilo,So^+ II^'at+I ||ilo,5o ^ 
+ |p(>^)||«o,5oi + £||Ap(^^)||«o,5o + ||AA(>^)||«„,5o + l|ASAr+i||«o,5o , 

which gives a bound on the norm of AZ. Notice that HCO-^^lj, ||if^>^^|| and || Ayl(>-^)|| are 
of order of e"^^°, namely of order once we set N such that N — [Kto/\ log A|]. We finally 
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define the constants Ci, C2, C3, C4 such that 

+ ||C^iV+l||«o,5o + ll^(>^)| 



Rq,So 


+ 


1 Gtv+i 


Ro,So 


< 




Ro,So 


< 


C2 


Ro,So 


< 




Ro,So 


< 





(52) 

With this setting we obtain: 

^||AZ||2 < {S + mr)\\AZ\\ + CiX^t + C2efxt + C3e + C4X^ . 

Based on the above formula and on || AF|| < || AZ||, we can draw the following conclusions: 

1. for some p\ > 0, || Ay || < pi for t of the order of if 

C2 = 0, namely if |bx ^^IkcSolk^-'^^lliio.So = 0' i-^- either |bx^^lko,So = or 

2. for some p2 > 0, ||AF|| < p2 for t of the order of the minimum between A~^ and 

(ep)-' if C2 ^ 0, i.e. \\p^^''^\\Ro,So\\s^^''mRo,So ^ 0. 

The two cases correspond to items %), %%) of the statement of the Theorem. More precisely, let 
us start with the case C2 — 0, i.e. p^^^ = or s^-^) = 0. Assuming that 5 + mr < aMv, 
C^e + C/^^ < iSMr'^, t < ^r^7Mr^ with Ci > 0, for some positive constants a, /3, 7, under 
the assumption that || AZ|| < r for some r > 0, we obtain 

^||AZ||2 < (a + /3 + 7)Mr^ , 

namely 

||Ay|| < ||AZ|| < V2(Q; + /3 + 7)r = pi (53) 
with pi < Rq. Taking into account (48) and (15), one obtains item %) of the Theorem, namely 

-1/(0)11 < 2CpA + pi fort < Ti = Coe^^° (54) 

with Co = (7Mr2)/Ci. 

Concerning item ii), since C2 7^ let cr > be such that for i < min(^^7Mr^, aMr"^), 
one has 

^||AZ||2 < (a + /3 + 7 + a)Mr^ , 
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namely 



I Ar II < II AZ|| < y/2{a + /3 + 7 + (7)r = p2, (55) 
with p2< Rq- According to (48) and (15), we obtain that 

||y(t)-y(0)|| <2CpA + p2 fort <T2 = min(^, ^) (56) 

with Co = i^Mr^)/Ci, C'o = {aMr^)/C2. 

When the M-convexity condition is violated (i.e. the second condition in (2)), by the assump- 
tion of quasi-convexity the first inequality in (2) must hold. Let T be either Ti or T2 as in (54), 
(56) with ||a;e||R(j replaced by 

C^= sup ||u;e(y)|| 

\\y-yo\\<P 

with p being pi or p2 as in (53), (55). Repeating the same argument as in [22], we assume that 
there exists an escape time Te such that || AZ|| — p and we show that this implies the inequality 

|a;e(y(s)) • AZ| < L||AZ|| VO<s<l. 

Then, using the same argument as for the convex case, we conclude that ||AZ|| < p, thus 
providing a contradiction. As before we have: 

\uje{Y{s)) ■ PAZ\ < {6 + mr)\\AZ\\ 
\ue{Y{s))-{Id.-P)AZ\ < CiX^T. 

Then, we have 

\uje{Y{s)) ■ AZ\ <{5 + mr) \\AZ\\ + CiA^ T . 
If T = Ti, p = pi, we obtain 

\uje{Y{s)) ■ AZ\ < ((5 + mr) II AZ II +CiA^ri 
< aMrpi + 'jMr^ 



V2(a + /3 + 7) 2(a + ^ + 7) 
< LWAZW , 



if 



+ 2(^;^)M 



^ y'2(a+/3+7) 2(a+/3+7) > 

If T = T2, p = p2, assume that T = [aMr'^) / {C2£p) (otherwise we recover the case T = Ti). 
Then, 

\u:e{Y{s)) ■ AZ\ <{5 + mr) \\AZ\\ + CiA^ . 
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Being 
we obtain: 



< 



\uJe{Y{s)) ■ AZ\ < {S + mr)\\AZ\\ + - aMr^ 

(7 

< aMr II AZ|| +7Mr^ 



V2(« + /3 + 7 + ^) 2(a + /3 + 7 + a)^ 



which is satisfied if the following condition holds: 

L 



Pi < 



^ ^2(a+/3+7+(T) 2(a+/3+7+(T) - 
□ 

Remark 5. Since we do not claim the result for any yo e A, but only locally under the 
conditions (4) and (5), we do not need to cover the whole phase space and therefore we do not 
need the analysis of the geography of the resonances, as it is usually done (see, e.g., [22]). 



4 Applications of the normal forms 

As we have seen in the Theorem, the stability time depends on the expressions of the terms 
p^x^^ s*^-^) appearing in the normal form equations, which represent, respectively, the 
conservative resonant part of the action variables and the dissipative resonant part pertaining 
to the angles, including the contribution of the modified frequency. In this Section we analyze 
several different examples, which well represent all possible situations which can be obtained 
with different choices of p^x^^ ^^'^ s^-^\ We illustrate these models with a twofold goal: to 
provide examples of cases i) and ii) of the Theorem and to illustrate an explicit evaluation of 
the resonant normal form. Since we do not aim to obtain stability estimates, we limit ourselves 
to the computation of the normal form in the non-extended phase space, i.e. in the variables 
X and y only. The experiments performed in this Section will be validated by the theoretical 
results of Section 5, where the estimates of the Theorem will be applied, showing linear as 
well as exponential stability times. 

All examples considered in the forthcoming Sections 4.1^.4 will have the following simple 
form: 

i = y + iJ-foiix^t) 
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y = -ehio,^{x,t) - fj.{y -rj) , 



where /oi and hio are periodic functions. In this case it is easy to decide which of the conditions 
i) or ii) of the Theorem are satisfied. Since we shall not need to consider Fourier modes less or 
equal, or greater than K, we drop the superscript by writing px,PY, sin place of Px^\ 
5(<^) Then, we can state that px = whenever the resonant part of hi^^^ is zero, otherwise 
px is different from zero. Concerning the function s, we can state that if the resonant part of 
/oi is not zero as well as if products of the form (/oi)™(/iio,a;)" with Q < m + n < N generate 
resonant terms or zero average terms of order /x, then the function s is different from zero. 
If the products (/oi)"^(/iio,a;)" with 0<m + n<A'^do not generate resonant terms or zero 
average terms of order ji, then s = up to the order N. 

4.1 Linear stability: case 7^ 0, s 7^ 

We consider the one-dimensional, time-dependent vector field given by 

X — J/ — ^(sin(,T — t) + sin(,T)) 

ij — —e{sin{x — t)+sm{x))—ii{y — ri). (57) 

Following the calculations of the proof of the Resonant Normal Form Lemma, the conservative 
transformation up to second order is defined by 



iJio{y,x,t) = 



sin(a;) 

y 

sin(2x — t) sm{t) sin(2x) 



2y2(2y-l) 2y2 gy^ ' 
while the dissipative contribution is given by 

(3oi{Y,X,t) = 
aoi[Y,X,t) = — 

cos(2X ~t) _ {2Y + 1) cos(i) cos(2X) 
au[y,X,t) - -2y2(2r -1) W'^ + 8^3 

MY,X,t) = 

ao2{Y,X,t) = + ■ 
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By choosing 

the normal form equations become 

X ^ Y 
Y - 



v{Y)^Y+ — + Os{e,i,), 



2^3 2Y 
—e sin(X — t) + Oz{e, n) 



where we recognize that px{Y,X,t) = sm(X - t), py(r,X,t) = -e/{2Y^), s{Y,X,t) = 
— sm(X — t) — ^. The Hamiltonian function in the extended phase space with U conjugated 
to time, associated to the normalized equations for = 0, is given by 



niY,x,u,t) 



Y2 ^2 

T + ^ + iy^ 



— e cos(X — t) + Os{e, ji) 



Replacing the normalized equations into the total derivative of %, one gets 



dn{Y,x,u,t) 

dt 



1 



1 



(^1 - cos(2X - 2t)^ + 03{6,fx) . 



(58) 



A typical orbit is shown in Figure 1 , where we integrate the normal form equations for e = 
IQ-^ and /i = 10"^ ^^jh initial conditions X(0) = and Y{0) = 1 + 6^. The left panel 
of Figure 1 shows the lift of {X, Y) to the universal coverage, while the middle panel shows 
the orbit back-transformed to the old variables {x,y). The dynamics starts on a rotational 
regime and drifts downwards; then it spirals along librational invariant curves until reaching the 
attractor. The right panel of Figure 1 provides the variation of the derivative of the normal form 
Hamiltonian, which tends to zero as the orbit reaches the attractor. The behavior is justified by 
(58) as the resonance is approached. 




10000 20000 30000 40000 50000 60000 



2 3 4 5 6 



2000 4000 



8000 10000 



Figure 1: Case px 7^ 0, s ^ associated to (57) for e = lO^^^ ^, = iq-3 ^j^^j fQj- jj^g initial 
conditions X(0) = 0, F(0) = 1 + Qy/e. Left: the lift of the normal form variables {X, Y) to the 
universal coverage. Middle: the trajectory in the original variables {x, y). Right: the variation 
of the derivative of the normalized Hamiltonian. 
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4.2 Linear stability at higher orders: case 7^ 0, s 7^ 

We consider the vector field 

X = y — fism{x) 

y — —e{sm{x — t) + sin(a;)) — ii{y — rj) . (59) 
The conservative normal form is defined by 

sin(x) 



y 

sin(25 — t) sin(t) sin(2a;) 



2r{2y - 1) 2p 8y3 
The dissipative transformation becomes: 

cos(X) 



aoi{Y,X,t) 



Y 



R (YXA - M^^-t) . sm(t) sin(2X) 

^ {l-3Y)cos{2X-t) (l-y)cos(t) cos(2X) 

' ' ^ 2y2(2y-i)2 ^ ^ gy^ 

Po2{Y,X,t) = 

, sin(2X) 

<^02{Y,X,t) = 

Higher normal form terms associated to (59) can be obtained in a similar way. For this model 
resonant terms occur at higher orders; for this reason we provide the following third order 
normal form equations: 

^ _ Y-^-I^ £^(2 - 5Y) cos(X - t) {6Y^ + 2y - 1) sin(X - t) 

~2Y^~2Y^ 2(i-2y)2y5 2(i-2y)2y3 
1 /I 2y \ 

+ 2^^' [y^ ~ {1 - 2Yy ) """'^^ ~ + 

. / X £^ sin(X — t) ^ , 
Y = -ssm{X-t)+ gy5_4y4 

The rj is the same as in the previous Section, the Hamiltonian associated to the normal form in 
the extended phase space is given by 

^, N Y'^ X e^cos(X-t) ^ , 

n{Y,X,U,t) = — + U+^-s cos(X -t)+ ^y4(2V-l) + ' 
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while the derivative of the Hamiltonian becomes 

dn(Y, X, U, t) _ eii^ sm{X - t) 



dt 



2Y 



The nonnal form produces a resonant term at third order. As a consequence, we observe a drift 
of the action variables, but on longer time scales. 





2000 4000 6000 8000 10000 



Figure 2: Case 7^ 0, s 7^ associated to (59) for e — 10~^, ji — 10~^ and for the initial 
conditions Xifi) = 0, 1^(0) = l + 6y/€. Left: the lift of the normal form variables (X, Y) to the 
universal coverage. Middle: the trajectory in the original variables {x, y). Right: the variation 
of the derivative of the normalized Hamiltonian. 



4.3 Exponential stability: case = 0, s 7^ 

We consider an example for which the normal form equations provide px — 0, but s 7^ 0. To 
this end, we modify the conservative part, so that the actions do not contain a resonant term at 
first order: 

X — y — /i(sin(x — t) + sin(x)) 

y = — e(sin(x — 6t) + sin(x)) — iiiy — r]) . (60) 
The conservative transformation to second order is given by 

, _ ^ sin(x — 6t) sin(x) 
y-6 y 

sin(2x) sin(2x — 6t) sin(2x — 12t) sin(6t) 



il^2o(y,x,t) = - ^ . . , ^^^^ - .^^ .„ + 



4y3-36y2 + 72y 8{y-6f 72y - 12y^ 
The dissipative transformation to second order takes the form 

/3oi{Y,X,t) = 

N cos(X) 
aoi{Y,X,t) = ^ 
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/3n(F,X,t) 



sin(2X - 7t) sm(2X - 6t) sin(2X) 



4y2 _ 38F + 84 4F2 _ sqy + 72 4F2 
^ sm{2X -t) ^ sin(t) ^ sin(5t) ^ sin(6t) 



2y-4y2 2y 10(F-6) 12(F - 6) 

, cos(2X) cos(2X - 7t) cos(2X - m) 

au{Y,X,t) = Q.^. + 7T7^ — ^T^N^TiT — ^ + 



+ 



8y3 2(7-2y)2(y-6) 8(y-6)(y-3)2 

cos(2A' -t) (y + 2) cos(f ) (y - 16) c()s(5t) (y - 18) cos(6t) 



2(i-2y)2y 2y2 50(y- 6)2 72(y-6)2 

MY,X,t) = 
The resulting normal form up to second order becomes 

2 2 2 

X = y-— ^ — --^--^-iism(X-t) + 03(e,ii) 

2(y-6)^ 2y'^ 2y ^ ^ ^ ^ '^^ 
y = 03(5,^), 

whereas the drift function is given by 

The Hamiltonian function in normalized variables corresponding to = in the extended 
phase space turns out to be 

y2 fy2 - ey + 18) 
n{Y,x,u,t) = - + T+ ^ 2(y-6)2y2 +Q3(^>^)' 

the time derivative of the Hamiltonian under the dissipative flow becomes 

dn{Y,x,u,t) 

It = ^3(e,/^), 

which shows the preservation of the energy up to the third order. Figure 3 displays the behavior 
of the lift of {X, y ) to the universal covering, the plot in the original variables and the graph of 
the derivative of the Hamiltonian versus time. The result shows that the dynamics takes place 
on an adiabatic quasi-periodic solution, which is consistent with the theoretical expectation. 

4.4 Exponential stability: case px 0, s = 

As an example which generates a normal form withpx 7^ 0, s = 0, we consider the differential 
equations 

X — y — /isin(6t) 

y — — £(sin(a; — t) + sin(a;)) — ii{y — rf) . (61) 
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500 1000 1500 

X 

Figure 3: Case px = 0, s ^ associated to (60) for e = 10^'^, /i = 10^'^ and for the initial 
conditions X(0) = 0, Y(0) = l + 6y/e. Left: the lift of the normal form variables (X, Y) to the 
universal coverage. Middle: the trajectory in the original variables (x, y). Right: the variation 
of the derivative of the normalized Hamiltonian. 
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The conservative transformation is given by 

/ /- ~ \ sinfx) 
Wio{y:X,t) = — — 

y 

while the dissipative transformation takes the form 

(3oi{Y,X,t) = 
aoi{Y,X,t) = --cos(6i) 



6 

sm{X - Qt) sm{X - It) siii(A' + 5t) sin(X + Qt) 
12Y - 2Y^ ~ 84 - 12Y ^ 12Y + 60 ~ 2Y{Y + 6) 
_ (3 - y) cos(X - 6t) (y + 3) cos(X + 6t) _ cos(X - 7t) cos{X + 5t) 

"in ' - (y - 6)2y2 + y2(y + 6)2 i2(y-7)2 " i2(y + 5)2 

l3o2{Y,X,t) = 
ao2{Y,X,t) = 0. 

The normal form equations are given by 



X = y-^ + 03(£,//) 

y = -£sin(X-t) + 03(e,/x) , 

with the drift function provided by rjiY, X,t) = Y — -^sii sin(X — t) + 0^{e, ji) . Note that we 
produce linear conservative resonant terms in the actions, but no resonant dissipative terms in 
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the angles. The Hamiltonian function associated to the normal form equations in the extended 
phase space becomes 

n{Y,X,U,t) = ^ + T + ^-s cos(X -t)+ Os{s, /.) , 
while the time derivative of the Hamiltonian flow becomes; 



dn{Y, X, u, t) 
Jt 



yielding the preservation of the Hamiltonian up to the normalization order. Figure 4 shows the 
behavior of the lift of {X, Y) to the universal covering, the graph in the original variables and 
the plot of the derivative of the Hamiltonian versus time. Also in this case, the result shows 
that the dynamics takes place on an adiabatic quasi-periodic solution, which is consistent with 
the theoretical expectation. 



iiiiiiiiiniiiiii 



nniiiiii 
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Figure 4: Case 7^ 0, s = associated to (61) for e = 10"'^, fi = 10"'^ and for the initial 
conditions X(0) = 0, Y{0) = l + 6^/s. Left: the Uft of the normal form variables {X, Y) to the 
universal coverage. Middle: the trajectory in the original variables (x, y). Right: the variation 
of the derivative of the normalized Hamiltonian. 



5 Application of the stability estimates 

In this Section we implement the Theorem to obtain estimates on the variation of the actions as 
given in Section 3. Let us fix the initial data as well as tq, sq (and related domain's parameters), 
K, 5. We assume that the frequency satisfies (4), (5) with a determined by (5). The smallness 
conditions on the parameters e, ji, say £ < £0, /j^ < jJ^o, come from (20), (23), (25), (28), (30), 
(33), (38), (39), (41), (43). 

We define the constants Cq and Cq as in (46), (47) and we set Cy = Cg + Cq- We recall that 
To, N, K are related by the expression 

To = — I logA| . 
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-1 
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-3 
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-3 
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10"=^ 


2.158-10- 


-4 


C2 


5.056 ■ 10" 


-3 


3.059 


• 10- 


-5 












C3 


2.01 




2.01 






3.208 - 


10-5 


2.01 




C4 


3.292 • 10- 


-5 


3.292 


• 10- 


-5 


1.006 - 


10-5. 


3.283 - 10- 


-6 


IIAFII 


2.408 • 10" 


-2 


2.408 


• 10- 


-2 


9.369 - 


10-3 


2.408 - 10- 


-2 


l|Ay|| 


2.421 ■ 10" 


2 


2.421 


■ 10- 


-2 


9.521 - 


10-3 


2.421 ■ 10- 


-2 


T 


2.692 • 10^ 




4.43- 


10^ 




1.699 - 


10^0 


3.309 - 10^ 





Table 1 : The main quantities of the Theorem for the examples of Section 4 from the remainder 

of a third order normal form. The parameters and initial conditions for all columns are: xq = 0, 
I/O = 1.01, ro = 0.05, fo = 4.9-10-2,f^, = 2.45- 10-^, i?o = 2.4-10-^, sq = 0.1, Sq = 5-10-^, 
^0 = 2.5 - 10-3, K=2Q,5 = 0.01. 

From (15) we determine Cp, while the constants Ci, C2, C3, C4 are computed as in (52). 
Table 1 provides the main quantities involved in the Theorem through the application of a third 

order normal form in the extended phase space. In particular, it provides the variation || AF|| 
of the normalized variables, the variation ||Ai/|| = \\y{t) — |/(0)|| of the original variables and 
the stability time T, which perfectly agrees with the theoretical result (linear or exponential 
stability time) of the Theorem. 

The results have been validated by a numerical integration of the equations of motion. Due to 
computer limitations, for the cases described in Sections 4.3 and 4.4 we had to stop to a time 
at most equal to 10^. Up to such integration times the numerical results are in full agreement 
with the analytical results. 
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6 Appendix A 



We briefly review the conditions which must be satisfied by the parameters e, n, so that the 
transformation from original to intermediate variables, as well as that from intermediate to 
final variables can be inverted; moreover, we provide conditions on the parameters so that the 
non-resonance conditions in the intermediate and final variables are satisfied. Compare also 
with [8] and [12]. 

6.1 Inversion of the conservative transformation 

With reference to (31), we invert the first transformation as 

x = x + r(^'^)(y,x,i) (62) 

provided that 

70||V'ri|fo,.oe'^°'^o-'<l, 

with 

II J- llr-o.so — II rj/ lko,so ■ 

for fo < ro, 5q < sq, sq = sq - Sq. 

6.2 Non-resonance condition after the conservative transformation 

Taking into account (5), we want that the non-resonance condition is satisfied in the interme- 
diate variables, say for a > 0: 

\uj{y) ■ k + m\ > ^ , \k\ + \m\<K. (63) 

The second of (31) can be inverted as 

y = y^eR^''\y,x,t), (64) 
for a suitable function R^^^ provided 

70||V'f^lk.o:^<l, 

for Tq < fo with 
Then we have 

yjjyy) • /c + m > - , 



2 



if 

a 

e < 



2K\\RW\\,,,, 
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6.3 Inversion of the dissipative transformation 

With reference to (12), the first equation can be inverted provided 

70||a(^)||,„,,„ e2^'°4-i<l, 
where Sq < sq. Inverting the equation as 

x^X + A(^'''\y,X,t), 

we have 

11^ llfo,so-6o — 11" lko,so • 

Thus we invert the second of (12) as 

y^Y + A(y'''\Y,X,t), 

provided 

To — rio 

with So < §0 — So, Ro < ro, being 

II A(^'^)||^„5o<ll>l(^'^^ 11.0,50- 

Notice that A^^'^^ can be bounded as 

11^^^'"^^ 11.-0,50 < \mo,o + ll/3f ^lko..-oll^^^'^^ll.-o.5o ■ 

Similar for the third equation in (12). 

6.4 Non-resonance condition after the dissipative transformation 

We now turn to the fulfillment of the non-resonant condition in the new set of variables 

\uj{Y)-k + m\>0, \k\ + \m\<K. 
Through the transformation 

and using (63) one finds 

\oj{Y) • A; + m| > ^ , 

provided that 

i^KILII/3^^^lk.o<^- 
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7 Appendix B 

From properties of analytic functions one can prove the following result (see also [8]) on the 
decay of the tail of the Fourier series. 

Lemma B.l. Let f — f{y, x, t) he an analytic function on the domain Cr^iA) x Cs(,(T^+^). 
Let f>^{y,x,t) = ^Q,rn)ez(+\\j\+\m\>K fjmiy) e^(^-^+"^*) and let < ctq < sq. Then, there 
exists a constant Ca = C'a((7o, K), such that 

WF^'lUso < a||/||.„,„+.,e-(^+^)-° , (65) 

with 

C..e(-«)*(i±^)'^'. (66) 
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